Canonical transformations of the extended phase space, 
Toda lattices and Stackel family of integrable systems. 



A.V. Tsiganov 

Department of Mathematical and Computational Physics, Institute of Physics, 
St. Petersburg University, 198 904, St. Petersburg, Russia, 
e-mail: tsiganov@mph.phys.spbu.ru 

Abstract 

We consider compositions of the transformations of the time variable and canonical 
transformations of the other coordinates, which map completely integrable system into 
other completely integrable system. Change of the time gives rise to transformations 
of the integrals of motion and the Lax pairs, transformations of the corresponding 
spectral curves and i?-matrices. As an example, we consider canonical transformations 
of the extended phase space for the Toda lattices and the Stackel systems. 
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1 Introduction 



It is well known, in classical mechanics any canonical transformation of variables maps a 
given integrable system into other integrable system. In this paper we consider compositions 
of the change of the time variable t and canonical transformations of the other coordinates, 
which map completely integrable system into other completely integrable system. 

On the 2n-dimensional symplectic manifold A4 (phase space) with coordinates {pj, qj}^=i 
let us consider hamiltonian system determined by the Hamilton function H{p, q). By defini- 
tion canonical transformation of the phase space Ai have to preserve canonical form of the 
Hamilton equations 

• _ 9H{p, q) ^ _ dH{p,q) 

Hi o ■> Pi o • l-*---*-/ 

OPi oqi 

As sequence, the action integral 

^ ~ j ^5-^ — H dt 

differential form YL^=i Pi dqi on M. and the Poisson brackets on M. are invariant with respect 
to the canonical transformations of the phase space M.. 

At the Hamilton equations (|1 . 1| ) and at the canonical transformations of M. the time t 
plays a role of parameter. If we want to consider change of the time t, we have to add new 



coordinate qn+i = t with the corresponding momenta Pn+i = H to the phase space A4 



The resulting 2n + 2-dimensional space A4 e in|, ^Sf is so-called extended phase space of the 
hamiltonian system. 

Canonical functional S on M.e has the following completely symmetric form 

^Pi Q-'i dr. 
1 i=i 

On the extended phase space M.e the Jacobi, Euler-Lagrange and Hamilton variational 
principles 5S = Q are differed by an additional constraint 

rL{pu. . . ,pn+i]qi ...,g„+i) = 0. (1.2) 



Here Ti is called generalised Hamilton function [|Il|]. According by (|1.2|) , any hamiltonian 
system is a conservative system on the extended phase space J^e- 

As an example, the Hamilton principle with generalised hamiltonian 

n = pn+i -E = H-E 

gives rise to the Hamilton equations on the extended phase space Ad e 

■ dn . dn . ^ 

Qi = T;-, Pi = -^, i = l,...,n + l. (1.3) 

opi dqi 



On zero- valued energy surface H = these equations are the initial Hamilton equations ([Ll] 
on the phase space M. and two additional equations 
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Below we shall consider conservative hamiltonian systems at dH/dt = only. 

By definition canonical transformations of the extended phase space A4 e preserve canon- 
ical form of the Hamilton equations (|1.3| ) or the Hamilton- Jacobi equation 

f + i/^O. (1.4) 

Invariance of these equations leads to the invariance of the canonical differential form 
Yl^i=iPidqi with respect to canonical transformations of the space M.e |]n|- Thus, change 
of the time variable 

t t , dt = v{p,q)dt (1.5) 
and transformation of the corresponding momenta H 

H ^ H=-^ (1.6) 

are canonical transformation of the extended phase space. Here H and H are variables of the 
extended phase space, but simultaneously they may be considered as the functions on the 
initial phase space [0. In general relativity function v{p, q) ( |1.5| ) is so-called lapse function, 
which determines transformation from physical time to coordinate time [l^ . 



Canonical transformations of the extended phase space M.e may be defined for any 
hamiltonian system. However, among all the hamiltonian systems a set of the completely 
integrable hamiltonian systems attracts a special attention. It is known, this set is invariant 
by canonical transformations of the initial phase space M.. In this paper, we want to discuss 
the following 

Problem: How to construct canonical transformation of the extended phase space A4e, 
which maps a given integrable system into the other integrable system. 

An exact technical definition of a completely integrable system is provided by the well- 
known Liouville theorem. Briefly, such system possesses n functionally independent integrals 
{Ij}^^^ of motion in the involution 

{Ij,Ik} = 0, j,k = l...,n, H = Ii. 

Thus, to construct canonical transformation of A^^; preserving integrability, we have to 
explicitly determine new algebra Ai of integrals by using a given algebra Ai of integrals of 
motion. 

For an integrable system we can consider dynamics with respect to the time tj conjugated 
to the integral Ij. Thus, in generic, we can add to the phase space M. all the integrals of 
motion {Ij} with the conjugated times {tj}. The geometry of the 4?T,-dimensional extended 
phase space Mlj = J^Q)j {Ij, tj} may be used to study of the transformations of the integrals 
of motion and conjugated to them times tj. 

By the Liouville theorem we can introduce the action- angle variables and 
in the vicinity of the common level surface of integrals 

Ma = {z eM : Ij{z) =aj, j = l,...,n.} 

Transformation of the time t ( p.. 51 ) acts on the half of the equations of motion only 

d<^j _ dipj _ Uj{si,...,Sn) 
—— — UJj[Si, . . . ,Sn ^ — ~- — 7- 7 , 

dt dt v{t,s,(p) 
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whereas other equations are invariant 



dt 



dsj 
dt 



0. 



If the lapse function vi p{s, q(s, if) ) depends on the action variables only 



v{s,ip) = v{s), 



change of the time (|1.5| ) preserves canonical form of the equations of motion in the action- 
angle variables 

5i = o. = (1.7) 

dt dt t;(si,...,s„) 

Moreover, at f = v{s) canonical transformations ( |1.5| - |1.6| ) of the extended phase space A4e 



map a given integrable system into other integrable systems. In this case mapping (|1.6| ) 
determines factorization of the initial Hamiltonian 



H{s) = v{s)-H{s) 



[1.8) 



up to canonical transformations of the action-angles variables {s,(p}. 

Thus, by the Liouville theorem we could construct canonical transformations of the 
extended phase space, which map a given integrable system into other integrable systems. 
However, it may of course be quite difficult to construct the action-angles variables for a 
given mechanical system, even if it is known to be completely integrable. Below we shall 
construct canonical transformation ( |1.5[ - |1.6| ) in the physical variables only. 

Change of the time preserving integrability induces transformation of all the machinery 
developed for integrable systems. Let us recall that the key idea which has started the 
modern age in the study of classical integrable systems is to bring them into Lax form 

{H,L} = [L,M], 

see reviews |]TB|, |TB|. The first matrix L or, more precisely, the coefficients of its characteristic 
polynomial -P(A) = det(L — A) determine integrals of motion. Change of the time gives 
rise to an algebraic transformation of the Lax matrices, a geometric transformation of the 
corresponding algebraic curves and transformations of the separation of variables methods. 
By using all the possible algebraic and geometric transformations, we can try to construct 
canonical transformations of the extended phase space preserving integrability. 

Transformations of the time variable and extensions of the phase space have been dis- 
cussed many times in classical and quantum mechanics. The wittingly incomplete list in- 
cludes: 

- solution of the equations of motion, as example for the Kowalewski top 

- the Kowalewski-Painleve analysis [|1^; 

- study of singularities of solutions equations of motion |T1]; 



- geometrical theory of integrable system [|Tl|, ^ 

- qualitative theory of dynamical systems H; 

- construction of the bi-Hamiltonian M and quasi bi-Hamiltonian systems 



the Birman-Schwinger formalism in quantum mechanics [17 
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As a first example, we consider integrable system determined by a natural Hamiltonian 

H{p,q) = Ho + V = ^g^jPiPj + V{qi, ...,qn), (1-9) 

where Ho{p) and V{q) are kinetic and potential parts of the Hamilton function. For any 
given energy E of the system we can use the Hamiltonian Ti = H — E to represent dynamics, 



provided that we impose the constraint H = 0, |23[. The passage to the geometric 



representation is accomplished by the Jacobi time transformation at 

v{q) = E- V{q) 

such that _ 

dt 

dt 



E - V{q) ' 

which maps orbits of an energy surface H = into geodesies of the Jacobi geometry, i.e. 
into the orbits of the following Hamiltonian 

H = ^ _ y^^-^ = J2 PiPj ' 9ij = {E- V)gij . ( 1. 10) 

The new Hamiltonian H will then give the same equations of motion on the surface 7i = 1. 
Here for a positive-defined kinetic energy Hq the term {E — V) is always non-negative in the 
physically allowed region. A modern discussion of this canonical time transformation may 
found in . 



Of course, the geodesic motion may be replaced by other systems. For instance, let us 
consider two-dimensional oscillator 



H = Ho + aV{q) = {pl+pl + b) + a {ql + qj] 



In this case free Hamiltonian Ho{p, q) does not pure kinetic part of H (|1.9| ). This choice may 
be considered as a shift of the energy surface Ti. = —b. However, it is more convenient to 
consider the free Hamiltonian Hq as a sum of the kinetic energy and the constant potential 
U{q) = b. For this system the Kepler canonical transformation of the time variable ( |1.6| ) 
with the function 

v{p, q) = V{q) =ql + ql 

preserves integrability. This change of the time t has been known by Kepler, the corre- 
sponding transformation of the Hamilton function has been studied by Levi-Civita |]T2[ and 



extended in [14 



After change of the time ( |l.(j| ) and the point canonical transformation of the other vari- 
ables (pi, gi,P2, 92) {p x,x,Py,y) the orbits of the oscillator maps into the orbits of the 
Kepler problem 

S- Ho{p,q) 2,2, b 

H = -— + a = p^+p+^== + a. 1.11 

v{q) 

Both these systems are degenerate Stackel systems in plane. We can construct at least four 



different Lax matrices for them. According to ^^e Kepler change of the time induces 
the following transformation of the Lax matrices 

L(A) ^ Z(A) = L{X) +(~ ° ] . (1.12) 
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The Lax matrix L{X) associated with the oscillator is defined on the initial phase space Ai. 
The Lax matrix L{X) associated with the Kepler problem is defined on the extended phase 
space Me- This transformation may be considered as a shift of the initial Lax matrix L{X) 
on the matrix from the extended phase space Me- Additional term in ( p.. 121 ) is the constant 
matrix with respect to the new time. 

Here both Hamiltonians H and H linearly depend on parameters a and b. Let the 
origin integrable Hamiltonian H{p,q; ai, . . . ,ak) linearly depends on arbitrary parameters 
(ai, . . . , afc). Then the function v{p, q) (|1.6|) may be founded by using additional propositions, 
for instance, that other integrals of motion are polynomials on parameters In this 
case the passage to the new time may be considered as a coupling constant metamorphosis 
acting on the one of the constant The similar additional propositions Kepler used by 

investigations of the celestial mechanics lows . 

It is known, in quantum mechanics, the eigenvalue problem of the Hamiltonian operator 

H 

H{p, q)^ = {Ho + aV + b)^ = , 
may be joined with the eigenvalue problem of the charge operator a 

H{p, q)^ = {Ho + {h-E) V-^)^ = -a^ . 

In quantum mechanics, such duality of the two eigenvalue problems has been used by Fok , 
Schrodinger [^] and many other. In the Birman-Schwinger formalism function v{q) is called 
a "sandwich" potential ll^. The canonical transformation of the time variable (|L6|) is an 



analog of this duality. Recall, in quantum mechanics zeroes of the potential v{p, q) = V{q) 
and signs of the parameters a and h determine compactness or non-compactness of the Sturm 
operator H in the corresponding space . In the classical mechanics zeroes of the potential 
^(P) Q) = signs of the parameters a and h determine behavior of the system with 

respect to the inversion of the time P]. 

Now let us consider an action of the canonical transformations of the extended phase 
space on the algebra of integrals of motion. Let the phase space of the initial system may 
be modelled on the dual space of certain Lie algebra g. It is known, the set of integrals {/, } 
gives rise to the subalgebra Ai in the corresponding universal enveloping algebra U{q*). 
A classical system is called integrable if the commutant Ai of the Hamiltonian H in f/(g*) 
contains an abelian subalgebra of the necessary rank. Of course, different abelian subalgebras 
in U{q*) determine different integrable systems on g*. However, the modern representation 
theory does not allows us to construct and classify all the possible abelian subalgebras in 

The search of the canonical transformations (|1.6|) may be reformulated as the construc- 
tion of the abelian subalgebra in U{q*) by using information about the known algebra of 
integrals Ai. For this purpose algebra Ai is too little, but algebra f/(g*) is too much. By 
using some additional propositions, we want to extend the first algebra (or restrict the second 
algebra) . 

Let the integrals of motions be either the fixed highest order polynomials in momenta, 
either polynomials in parameters, such that 

Lk i,k 
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In the fist case instead of origin algebra Ai we can introduce new algebra A'^ D Ai generated 
by momenta and corresponding coefficients Jj'^ in ( |1.13| ). In the second case new algebra 

Aj D Ai is generated by parameters and coefficients Jj'^. Expansions ( 1.13 ) give rise to 
grading of the new algebras A^/'^ by power of momenta or parameters. 

By the Liouville theorem initial system is integrable, if the algebras Af"' contain the 
abelian subalgebras of the necessary rank. If one of the elements v{p,q) = or v{p,q) = 

P-^ is invertible in a given representation, we can introduce a completed algebra A^i"" = 
such that rank^j'" > rank^/. The search of the two abelian subalgebras in the 
algebra ^j'"} with elements H and H ( |1.6|) is equivalent to construction of the canonical 
transformation (|1.6|) . 

Next we can consider two integrable systems on a common phase space M.. By using 
two algebras of integrals Ai and Aj we can construct new algebra Aij = Ai ® Aj or 
Ajj = Ai ® Aj. Of course, algebra Ajj should be much richer than Aj and Aj, separately. 
This algebra has two abelian subalgebras, which associate with two initial integrable systems. 
If one of the elements f (p, q) = Ij or v{p, q) = Jj is invertible, we can introduce a completed 
algebra Au = {v~^,Aij} such that TankAu > rankAi. Now, existence of the third abelian 
subalgebra in Ajj with elements H ( |1.6| ) is provided existence of the canonical transformation 

(PD- 

Obviously, initial algebras Ai and Aj contain all the necessary information about the 
completed algebras ^f'" or Aij. Thus, we have to extract this information only. As an 
example of the solution of such algebraic problems, we consider integrable systems related 
by canonical transformations with the Toda lattices and the Stackel systems. 

In the second Section we propose canonical change of the time for the Toda lattices 
similar to the Kepler change of the time ( p..ll| ). Let us rewrite the Hamiltonian of the Toda 
lattice as a sum of the free Hamiltonian Hq and potential Vp 

where P is a system of simple roots and (3 is one of these roots. Here free Hamiltonian Hq is 
a sum of the kinetic part with the nontrivial potential. Similar to the Kepler transformation 
( |1.11| ), change of the time ( [1.5| , |1.6|) preserves integrability at 



v{q) = = ef^^"^ (3eP. (1.14) 



The corresponding mapping of the Lax matrices looks like ( |1.12|) . To construct these new 



integrable systems we shall use other automorphisms of the Ai subalgebras in the grading 
simple Lie algebra q associated to the Toda lattice. 

In the third Section we consider some pair of the Stackel systems with a common Stackel 
matrix, defined by two families of integrals {Ik}^=i and {Jk}^^i ■ In this case ratio of any 
two integrals H = Im/ Jm ( |1-6| ) and the inner product (I A J) of the vectors I and J on 
may be used to determine of the third integrable set of integrals of motion. To construct 
corresponding change of the time ( |1.6|) we explicitly used the Jacobi ideas in the separation 
of variables method. Recall, Jacobi proposed to construct various classes of completely 
integrable systems starting from a set of separated one-dimensional problems. 
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2 The Toda lattices. 



Before proceeding further, it is useful to recall some known facts about generalised Toda 
lattices (all details may be founded in review |l8l). 

Let g be a real, split, simple Lie algebra rankg = n. Let ) be its Killing form, let a be 
a split Cartan subalgebra, A the associated root system, A+ the set of positive roots and P 
the system of simple roots. For a G A+, let be the corresponding root space and G go, 
a root vector. It will be convenient to normalize in such a way that K{ea-, e_Q,) = 1. 

The root space decomposition, g = a + ®aQoL) gives rise a natural grading on g, the so- 
called principal grading. Let g+ = ©i>ogi be a Borel subalgebra of g and g__ be the opposite 
nilpotent subalgebra. The g = g+ + g_ is the generahsed Gauss decomposition for simple Lie 
algebra, which induces a dual decomposition g* = g+*+g+ and an injection i = g+ > g*. Let 
Oa C g+* denote the coajoint orbit of some a G g+*. The image i{Oa) G g* has, in general, 
nothing to do with the coajoint orbit (9i(a) of i{a) in g*. So i* of the Casimir function on 
g* can be a complicated, non-constant function on Oa- The Adler-Kostant-Symes theorem 
claims that any two such functions i*(f)i and 2*02 are in involution with respect to the natural 
symplectic structure on Oa- 

The Killing form K allows us to identify g* with g_j, so that 

g+ = ©i<oSi , 0- = ®i>o9i ■ 

Choose a vector 

a = tta e_Q , G M 

and let Oa be the g_|_ orbit of a in a + g_i. The points of Oa have the form 

^ = CattaC-a, CI, C„ > 0. 

The orbit Oa does not really depends on the a^'s but only on their signs. 

It is convenient to introduce the variable g G a such that Cq = exp a{q). Fix a basis {hj} 
in a, and let {h'^} be the dual basis with respect to the Killing form K{hi, h'j) = 6ij. Set 

Qj = K{q, hj) , pj = K{p, h'j) , p, g G a . 

Then 

{Pi^Pj} = {Qi, Qj} = , {Pi, Qj} = Sij , 

i.e. the variables {pj,qj} are canonical. The orbit Oa is parametrized by the canonical 
variables as follows 

n n 

i = ^Pihi + ^aa- exp(^ qiK{a, h[) ^ ■ e_„, . 

i=l aeP i=l 

Replace the inclusion i = q\ "—>■ Q* hj its shifted version i + e, where e is an element of g* 
which annihilates [g±,g±]. 

For the Toda lattice the orbit Oa is an orbit of g+, what we are really interested in the 
orbits of the full algebra gn = 0+ ©0-- Let us translate Oa by adding to it a constant vector 
e = XlaGP^" which may be regarded as a one-point orbit of g_. The resulting orbit 

Oae = Oa + e (2.1) 
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is parametrized by the canonical variables as follows 



L = ^ Pihi + ^aa- exp K{a, q) ■ e_a + ^ . (2.2) 

i=l aeP aeP 

Let us consider the simplest Hamiltonian on Oae generated by the Killing form on g 

HiX) = ^KiX,X). (2.3) 
Its restriction to Oae (|2.2| ) is given by 

H{p,q) = ^K{p,p) + J2ciae"^''K (2.4) 

To keep in mind consequent canonical transformation of the time, we choose the second Lax 
matrix at the following special form 

^L=[L,L„], L_ = - ■ expK(a,g) ■ e_„ . (2.5) 

Any finite-dimensional linear representation of g gives rise to a matrix valued function L 
( p. 21) on Ai, the coefficients of its characteristic polynomial are integrals of motion. 



The behavior of the dynamical system with the Hamilton function (2^) depends crucially 
on the signs a^. The associated hamiltonian flow is complete if and only if aa > 0. In this 
case, its solution is reduced to the generalised Gauss decomposition in the corresponding Lie 
group. Systems with other signs of Oq, are of course also solvable, but their solutions blow 

up ig. 

Now let us construct canonical change of the time variable for these Toda lattices. Recall, 
in a shifted version of the Adler-Kostant-Symes scheme in order to get orbit Oae in 7V1 ~ 
0p = 0+ ® 0- translate orbit Oa living in by adding to it a constant vector e from the 
remaining part of Ai. Let us replace the phase space A4 on the extended phase space A4e- 
In this case we can also translate the same orbit Oa in g*^_ by adding to it a constant vector 
from the remaining part of the whole space A4e- As above, this vector has to be a character 
and a constant with respect to the new time. The third condition is that the initial invariant 
polynomial (|2.3|) have to generate the coupling constant 

K{L,L) = -b 



instead of the Hamiltonian H, as for the Kepler problem [27 



Proposition 1 For each simple root (3 & P and for any constant 6/3 G M the following 
canonical transformation of the extended phase space A4 e 

dt = e'^^'i^ ■ dt , 

(2.6) 

Hp = e-^(«) -(^H + b^ 
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maps the Toda lattice into the other integrable system. This canonical transformation induces 
the following transformation of the Lax matrices 



L-H, 



(2.7) 



such that 



{i7,L} = [L,L_]. (2.8) 

Here H , L and L_ are the Hamiltonian ( ^.4\ ) and the Lax matrices Sj^. 5| j for the corre- 
sponding Toda lattice. 



As for the Toda lattices, restrictions of the invariant function on the orbit 10 give rise 
to the new set of integrals of motion. 

There are n = rankg functionally independent invariant polynomials on g. Restricted to 
the orbit Oae they remain functionally independent and give rise to integral of motion for the 
Toda lattice. In (p.7|) we translate Oae by adding to it a "constant" vector proportional to the 
element of the universal enveloping algebra. All the invariant polynomials are invariant with 
respect to this transformation. Thus, we can construct n independent integrals of motion in 
the involution for the system with the Hamiltonian ( p.6| ). 

The number of the functional independent Hamilton functions Hf^, j3 E P depends on the 
symmetries of the associated root system. For the closed Toda lattices associated with the 
affine Lie algebras canonical time transformation has the similar form. The associated with 
the Lax matrices L ( p.2| ) and L ( |2.7| ) spectral curves depend on the choice of a representation 
of 0. Therefore, geometric transformations of the spectral curves will be below considered 
on the some examples only. 

Of course, canonical transformations ( |2.6|) induce transformations of all the machinery 
developed for the Toda lattices. As an example, for the Toda lattices associated to the 
classical infinite series of the root systems An, Bn, Cn and Dn we can introduce another 



2x2 Lax matrices jTOf. For brevity, here we restrict ourselves the An root systems only. 
In this case, the second Lax pair representation has the following form 



T = T,T2---Tn 



where 



A 



-Pj 







such that 



-e-* 



{H, Tj} — Tj Aj - Aj_i Tj , Aj - 

Change of the time (|T]^, P^ ) associated with the root [3 = Sj — Sj^i by 

v{q) = e"^^'') = exp(gj+i - qj) 
gives rise to the following transformation of the Lax matrices 



(2.9) 



T = TiT2---Tn 
the Lax equations 



T = Ti---T,^i 



H + b 




T, 



i+2 ' 



Tn, (2.10) 



{H, T} = [T, An] ^ {H, T} = [T, v{q) A 
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and the spectral data of the Lax matrices 



detT = nLi«A:, trT = A" + A"-ip + A"-2 - ^) +•■• . 

det f = (a,- - H) ■ rifc^, afc , tr f = A'^ + A"-i P + A"-^ + 6^ + . . . . 
All these transformations look like corresponding transformations for the Stackel systems 



27, 28 



Let us consider some heuristic arguments concerning the Lax matrix transformation 
(p.8|). Recall, constructions of the orbit Oae and the corresponding classical P-matrix are 
closely related to the principal grading of g. This grading defines all the possible embedding 
of the three-dimensional subalgebra Ai ~ sl{2) into g. 

Let {e_, e_|_, h} be generators of the Lie algebra sl{2) 

[h,e_]=e_, [h,e+] = -e+, [e_,e+]=2h, (2.11) 

and the element 

A = h2 + l(e_e+ + e+e_) (2.12) 

of the universal enveloping algebra be Laplace operator in SL{2). Let us consider infinite- 
dimensional irreducible representation W of the Lie algebra sl{2) in the linear space V such 
that 

W: {e_,e+,h} ^ {e,/,/i} e End(y). 

Let operator e be invertible in End(V) and v^(A) be arbitrary function on the value of the 
Casimir operator ( |2.12| ), then the mapping 

e_ e'_ = e_, h ^ h' = h, e+ ^ e'_^_ = e+ + el^ • v^(A) (2-13) 

is an outer automorphism of the space of infinite-dimensional representations of s/(2) in V 
||26| . These mappings shift the spectrum of the Laplace operator A (|2.12|) by the rule 

A ^ A' = A + v?(A). 

In particular by v^(A) = — (A + 6) one gets A — > A' = —b. So, instead of the spectrum of 
the Laplace operator A on the group SL{2) one gets the spectrum of the coupling constant 
h. 

Now we turn to the Toda lattices related to the Lie algebra 0, which contains subalgebras 
Ai associated to the roots (3 & P. The maps (|2.8|) may be closely related to the automorphism 
( p.l3| ). Instead of the restriction A ( p.l2|) of the Casimir operator H{X) ( p.3|) on Ai one 
have to substitute restriction H{p,q) ( p.4| ) of the same Casimir operator on the orbit Oa,e- 
Note, the outer automorphism ( p.l3|) non-trivial acts on the one nilpotent subalgebra of s/(2) 
only. Extension of this mapping on the orbit Oae non-trivial acts on the constant one-point 
orbit of 0_, too. 

Of course, the more justified consideration of the canonical transformation of the time 
variable requires a more advanced technique of representation theory. 
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2.1 Examples 

Let us consider periodic three-particles Toda lattices associated to the affine Lie algebra 
C{As). The corresponding Hamiltonian H ( p.4| ) reads as 

H=^{pI+pI + pI) + ai + as e'?^-^^ + 03 e'^'-"' , (2.14) 

and the Lax matrix L{X, fi) ( p.2|) has the following form 

X-Pi aie«-'?2 /i-i 
L(A,/i) =1 1 \-p2 aae''^-''^ 

/i ■ 03 e^^""! 1 A - p3 



The Lax equations (|2.5|) is linearized on the Jacobian of following algebraic hyperelliptic 
curve 

C: detL(A,/x) = 0, (aiaaas/x + -) + A^ + A^ P + A f — - ifj + = . (2.15) 



Here complete momenta of the system P = pi +P2+P3, Hamiltonian H ( p.l4|) and additional 



integrals of motion K are polynomials of the first, second and third order in momenta, 
respectively. 

The canonical transformation (|2.6|) maps the origin Hamiltonian ( p.l4| ) into the following 
Hamiltonian 

H = e"^^-^' ■ {H + b) , beR. (2.16) 
The associated Lax matrix L ( p.7|) is given by 

/ 0^ 

L = L~ I a^^H I . (2.17) 

\ 

Its characteristic polynomial is 

C: detZ(A,/i) = 0, (aiaafls/i + ^^~'^ ) + X^ + X^P + \{— + b]+K = 0. (2.18) 

ai /i \ 2 



Here, in comparison with (|2.15| ), the linear in A term is proportional to the constant b 



instead of the Hamilton function H. The new Hamiltonian H is related to the second 
spectral parameter now. For the first time similar transformations of the Lax matrices and 
of their characteristic polynomials have been observed in p8| . 

The corresponding 2x2 Lax matrix for the three-particle Toda lattice is equal to 



T = TiT2T3 



X-pi ai \ / A - p2 (32 \ / A - ps e^^ 



Change of the time maps this matrix into the following matrix 

X- P2 a2 e''^ 

^1 ■ I -e-^^ (1 - ai'ij) • ■ ' 
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such that 

detT = aiaaag, tiT = + X'^ P + X - + K , 

detf = (01-^)0203, trf = A^ + A^P + A ^ ^) 

If we put aj = 1, the Poisson bracket relations for the initial 2x2 Lax matrices are closed 
into the Sklyanin quadratic r-matrix algebra |10 



{r(A) , r(/i) } = [r(A - /i) , r(A)T(/x) ] , r(A - /x) " 



A — /i 

Here the standard notations are introduced: 

T(A)=T(A)®/, T(/i) = / ® r(^) , 



and n is the permutation operator of auxiliary spaces Change of the time (|2.16| ) maps 
the Sklyanin r-matrix algebra into the following poly-linear algebra 

{r(A) , Tifi) } = [ r(A - , r(A)r(/x) ] + [ ri2(A, /i) , T(A) ] + [ r2i(A, /i) , T(/x) ] , 
with dynamical r-matrix 

ri2(A,/i) = A3(A)®aT3(/i), a = 



1 




r2i{X,fi) = -nri2(/i, A) n. 



Here T3(A) and A3 (A) were defined in (p.9|). 

By using complete momenta P, we can reduce initial 6-dimensional phase space to the 
4-dimensional space. After this reduction and after the following point canonical transfor- 
mation of the coordinates 

qi = ^(1 + iVS) Inx + ^(1 - iVS) \ny , 



(2.19) 



1 i 1 i 

— 1 ^)xpx (1 1=] 

2^ Vs' ^ 2^ 



the new Hamiltonian H reads as 



H =pxPy + — + a2 y'' + a3 x'^ y'' + a, , (2.20) 
xy 
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where Zj are roots of the quadratic equation {z — Zi){z — Z2) = 2;^ + 3z + 3 = 0. For the first 
time the system defined by Hi (|2.20|) has been found in [Q. 

In conclusion we present integrable systems related to the two-particle Toda lattices 
associated to affine algebras x'^^ . After an appropriate point transformation of coordinates, 
all the Hamilton functions have a common form 

H = PxPv + — + ax^' + ex"*' y""^ + d, a,b,c,d eR 
xy 

where zi^2 and si,2 be the roots of the different quadratic equations. Below we show these 
equations only: 

^2 + 3^ + 3 = s2 + 3s + 3 = 

2^ + 4^ + 5 = + 4s + 5 = 
^2 + 2^ + 2 = s2 + 3s + 5/2 = 

^2 + 2^ + 2 = s2 + 2s + 2 = 
^2 + 2^ + 2 = (s + 2)2 = 

^2 + 2^ + 4 = s^ + 5s + 7 = 
^2 + 2^ + 4 = s^ + 3s + 3 = 
^2 + 3z + 7/3 = s2 + 3s + 3 = 

The corresponding second integrals of motion K are polynomials third, fourth and sixth 
order in momenta. Note, for the algebra Ag^^ all the three Hamiltonians Hp, (3 G P are 
equivalent. Two different Hamilton function ( p.6|) associated with the algebras b!^\ Cg^'* 
and -02^''. For the G^2^ algebra we have three different Hamiltonians ( |2.6D . 

3 The Stackel systems 

In this Section we propose an interesting extension of the integrable family of the Stackel 
systems [^. As an example, we discuss here a family of the two-dimensional integrable 
systems in detail. In the two limiting cases, the corresponding systems possess the following 
Hamilton functions 

k 

(3.1) 

k 

H = pl+pl + a^xy) k+i , 

where a and k are arbitrary parameters. At A; = 1 the first Hamiltonian coincides with the 
Hamiltonian of the Kepler problem. At A; = 2 the second integrable Hamiltonian has been 
found by Fokas and Lagerstrom It is known, both these systems are dual to the some 
Stackel systems (see review ^ and references therein). By using similar duality we shall 
prove integrability of the general systems ( |3.1| ). 

Let variables {pj, qj}^=i be coordinates in the phase space M^"' with the standard Poisson 
brackets 

{Pj,qk} = 5jk, j,k = l,...,n. 



4': 
U2 ■ 
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Let us consider two integrable hamiltonian systems on the common phase space M^". These 
systems are defined by the two sets of independent integrals of motion {Jj}"^^ and {Jj}^^^, 
in the involution 

{Ij,Ik} = and {Jj,Jk} = 0, j,k = l,...,n. 

By the Liouville theorem for given integrable systems we can introduce two family of the 
action-angle variables {sj, V^j}j=i and {sj, for the each of integrable systems. Integrals 

of motion depend on the action variables only 

Ik = Ik{Si, . . . , Sn) , Jk = Jk{Si, . . . , S„) . 



According to expansion ( |1.8| ), if two a'priori different systems of the action coordinates are 
related 

Sj = Sj{Si,.. .,Sn) (3.2) 

then the canonical transformations of the extended phase space at 

v{p,q) = v{s) = Ikip,q) , or v{p, q) = v{s) = Jk{p, q) 

map a given pair of integrable system into the third integrable system. 

Thus, we have to verify condition ( |3.2| ) and have to construct explicitly the third family 
of integrals of motion in the initial physical variables. For instance, by using inner product of 
the two independent vectors of integrals I and J in we can introduce some antisymmetric 
matrix /C = (1 J). Any column or row of this matrix defines a set of the n — 1 independent 
functions 

^ij ~ (I ® J)ij ~ J j IjJi^ h j ~ ■ ■ ■ y ^ ■ 

Under some restriction on the initial integrals {Ij} and {Jj} we could construct from them 
a third integrable system on the same phase space. 

Proposition 2 If all the differences of integrals of motion {Ij — Jj) with the common index 
j = 1, . . . ,n are in the involution 



|/j - J, ,4 - Jfej = 0, J, fc = l,...,n, (3.3) 

then the ratio of integrals 

= 

J, 

and n — 1 functions Kj, j ^ m 



Km = ^ (3.4) 



i^i = ^= ^'"'^% =K^Jj-I,, m^j = l,...,n (3.5) 

are integrals of motion for new integrable system on the same phase space. 

By definition all the new integrals Km and Kj are functionally independent. Thus, the proof 
is straightforward calculation of the following Poisson brackets 

{Km, K,} = ( J,} + {Jm, IA)=-^ {l^n - Jm, Ij " = . 
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m 



and 

{K,,Kk} = K^{J„Kk}-{I„Kk} = 

= Jk [k^ {Jj, Km} - {Ij, Km}^ - Km [{Jj, h} + {Ij, Jk] 

= Jk{{Kj,Km} + Km[lj-Jj,Ik-Jk}=0. J^k^ 

Thus, canonical transformation ( p.4[ - p.5|) of the extended phase space 

(I,J)-.K 

preserves the property of integrabihty. To apply this transformation we have to find two 
integrable systems satisfying condition (|3.3| ). Below we prove that the Stackel integrable 



systems may be considered as the main example of the systems satisfying condition 

Let us briefly recall some necessary facts about the Stackel systems The non- 

degenerate n X n Stackel matrix S, whose j column s^j depends only on qj 

detS^O, ^ = 0, j^m 

oqm 

defines functionally independent integrals of motion {Ik}^^i 

n g 

^k = Yl ^J'^ W + ^^j) ) ' ^^'^ = detS ' ^^'^^ 
i=i 

which are quadratic in momenta. Here C = [cik] denotes inverse matrix to S and Skj be 
cofactor of the element s^j. 

We can see that in practical circumstances the Stackel approach is not very useful be- 
cause it is usually unknown what canonical transformation have to be used to transform a 
Hamiltonian ( ^.61 ) to the natural form H = T{pi, . . . ,pn) + V{qi, . . . ,qn). This problem was 
partially solved for the Stackel systems with a common potential Uj = U, j = 1, . . . ,n only 

According to if the two Stackel matrices S and S be distinguished the m-th row 
only 

^kj = Sfcjf ) k ^ m , 

the corresponding Stackel systems with a common set of potentials Uj and with the Hamilton 
functions Im and Im are related by canonical change of the time 

v{q) 

where _ 

det S(gi, . . . 



detS(gi,...,g„) 

The canonical transformation ( |3.7| ) connects two Stackel systems with the different Stackel 
matrices and with the common set of potentials Uj. 
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Let us consider a pair of the Stackel systems with a common Stackel matrix S and with 
the different potentials. Namely, in addition to the system with integrals {Ik} ( |3.6| ), we 
introduce the second integrable system with the similar integrals of motion 

n 

Jk = Yl Cjfc (P? + Wj{qj)) , k=l,...,n. (3.8) 

i=i 

Here even one potential Uj (qj) has to be functionally independent on the corresponding 
potential Wj (qj). 

Proposition 3 Any two integrable systems defined by the same Stackel matrix S and by 
the functionally independent potentials Uj (qj) and Wj (qj) satisfy the necessary condition 
( pT^ j of the previous proposition. Thus, the ratio of the two Stackel integrable Hamiltonians 
defines new integrable system 

Km ^ ^ (^m,! J-rri) ; ~^ • (3-9) 

It is obvious, all the integrals and Jk are differed by the potential part 

2 

(4 -Jk)=Y^ [ Ujiqj) - Wjiqj) ] 
i=i 

depending on coordinates {qj} only. Thus, systems with a common Stackel matrix S satisfy 
condition ( |3.3| ). 

The Hamilton function (|3.9|) has the following form 

E,=lCjm [Pj + Wj{qj)\ 

This Hamiltonian if is a rational function in momenta, but next one can try to use canonical 
transformations to simplify it. In rare case, one obtains again a natural type Hamilton 
function as will be shown below. 

As for the usual Stackel system, the common level surface of the new integrals (|3.4| ) 

M^ = [ze M^" : Kj{z) =aj, j = l,...,n} 

is diffeomorphic to the real torus. Namely, substituting 

Vj (qj) = (1 - am)-\ Uj (g,) - Wj (g,) ) , 

into the definitions (|3.10|) and ( p.6| , p.8|) we obtain the following equations 

2 



(3.11) 



2 



j=l 
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After multiplication of these equations by the Stackel matrix one immediately gets 



Pi 



^ ^ k=l 



where S{qi, . . . , qn) is a reduced action function The corresponding Hamilton- Jacobi equation 
on M„ 



dS_ dS_ dS_ 



0, 



dS_dS_ 

dqj dqj 



E. 



admits the variable separation 5(gi, . . . ,qn) = J2^=i where 



Sjiqj 



\ 



^^PkSkjiqj) -Vjiqj) dqj. 



k=l 



Thus, coordinates qj (t, ai, 



, an) are determined from the equations 
Skj{\)d\ 



k 



n, 



where the constants aj and Pj are related by 

Thus, the solution of the problem is reduced to solving a sequence of one-dimensional 
problems, which is the essence of the method of separation of variables. Next, the integration 
problem for equation of motion is reduced to solution of the inverse Jacobi problem in 
framework of the algebraic geometry p7 |. 

Recall, the first canonical transformation (|3.7|) may be related to the ambiguity of the 
Abel map. So, it would be interesting to investigate the underlying algebro-geometric origin 
of the construction (|3.9|) . 

For the some Stackel systems with uniform potentials Uj = U, j = l,...,n we can 
construct the 2x2 Lax matrices lESl. Of course, transformation of the time (U-7i) induces 



transformation of the Lax matrices. As for the Toda lattices, these transformations may be 
reduced to the adding a constant with respect to the new time matrix ^ to the initial 
Lax matrix. Transformations of the Lax matrices by new mapping of the time ( |3.9| ) will be 
published elsewhere. 



4 Examples 



As we have mentioned before, the Hamiltonian H ( |3.10| ) have a rather unusual expression. 



However, in some cases suitable canonical transformations can reduce it to a sum of the 
kinetic energy and the potential energy. Note, necessary choice of such transformations is a 
generic problem for all the Stackel systems. Thus, in this Section we present several concrete 
two-dimensional systems only. 

Let us consider polar coordinate system on plane with the usual coordinates {pr, r) and 
0) instead of (^1,2, 91,2)5 respectively. We take the first system with the axially symmetric 
potential 



h=pl + ^ -a^r^'' + b, h = Pc^ , a,b,neR. (4.12) 
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The second system is associated to a free motion 

2 

^l=Pr + ^, J2=V^- (4.13) 

These systems belong to the Stackel family of integrable system with the following common 
Stackel matrix 

1 

-r-2 1 

In this case Si2 = C21 = and h — -h = 0, it allows us to use the second integrals in the 
non-Stackel form ( p.6|) . 

The ratio ( p.4|) of the Hamiltonians ([4.12| - ^?T3|) may be rewritten in the form (|3.1| ) by 
using the set of canonical transformations. Let us begin with the usual transformation of 
the curvilinear coordinates to the cartesian coordinates 



r = \/uv , = z arctan 



u — V 



u + V 



(4.14) 



UPu + VPy . 

Pr = , P<p = i [UPu-Vp^] 



Then, one permute coordinates and momenta [u pu) and {v ^ p^) such that new Hamil- 
tonian ( |3.10|) becomes polynomial in momenta 



i/ = ^ = -^^ + ^ + l. (4.15) 

Ji A uv Auv 



In conclusion, we have to use the point canonical transformation 

1 k 

Px=PuU^+^, X = (1 + 1/fc) U'^+l 

1 k 

Py=PvV''+^, y = {1 + 1/k) v'<:+^ , 



which converts the Hamiltonian (|4.15| ) into the following form 



k 

H = plp'l + a{xy) fc+i 

after multiplication on a suitable constant and rescaling parameters. 

Note, the Stackel matrix S and the set of potentials Uj{qj) are determined on the half 
of the phase space M^" and depend on coordinates qj only. Thus, we have some freedom 
related to the different canonical transformations of momenta 

{Pl,... Pn) (Pl,...,p„), 

(4.16) 

Pi- Pi = 2 , p . + p . = . 

oqi 



Here F{qi, . . . qn) is a generating function of the transformations ( |4.16| ) depending on coor- 



dinates {qj}, which are invariant with respect to transformation (|4.16| ). 
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If condition (|3.3|) is invariant under these canonical transformations, we can apply them 
( [4.1(j| ) to construct new integrable systems by the rule ( |3.4] ). Although no general procedure 
exists for this, one interesting case is known. 

As above, one takes system with the axially symmetric potential ( |4.12| ) and system 
associated with a free motion with integrals 



Ji 



New momenta {pr 



Pr = Pr 



relate with old ones {pr 
d f{r) cos(?T,0) 



dr 



n 



J2=P4>- (4.17) 
I, ) by canonical transformation ( 4.16|) 

V =P0 + /W sin(n0). (4.18) 



Here /(r) be any function on variable r and n be arbitrary parameters. 

Both these systems belong to the Stackel family of integrable system associated with 
a common Stackel matrix. In this case /2 — J2 7^ and condition ( p.3|) do not invariant 
by transformation ( |4.18| ). Let the second integrals be square root from the usual Stackel 
integrals (|3.6| ). If this form of the second integrals are used, the pair of the systems (|4.12| ) 
and ( |4.17| ) satisfies condition ( p.3[ ) by 



fir) 
fir) 



n 
r 



fir) 



cr 



c G 



Let us get over to the cartesian coordinates ( [4.14| ) and conjugated momenta 

UPu + VPy 



Pr 



I [UPu - Vp^ 



(4.19) 



After permutation coordinates and momenta {u <-> pu) and {v ^ p^) (analog of the Fourier 
transformation in the quantum mechanics), one gets 



H 



h 
Ji 



1 



(c'K-^r' - a'ptpt - 2 c {vp:-' + up:-' )) + 



(4.20) 



4-ut> 4Mf 
At c = we discuss this system before ( [4.15|) . Now, we consider the second limiting case by 

a = c, k = n — \ , 

that simplifies potential part of the Hamiltonian ( |4.20| ). 

As above, the point canonical transformation ( 4.16| ) converts the Hamiltonian ( 4.20| ) into 
the following form 

k 

H =p^,+p'; + a{xy) k+i 

after multiplying on a suitable constant and rescaling parameters. 

Thus, we present a family of two-dimensional integrable systems, which includes the 
Kepler and Fokas-Lagerstrom potentials simultaneously. 
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